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ABSTRACT 


We have solved for the potential flow downstream of the terminal shock of the solar wind in the limit of 
small departures from a spherical shock due to a latitudinal ram pressure variation in the supersonic solar 
wind. The solution connects anisotropic streamlines at the shock to uniform streamlines down the heliotail 
because we use a non-slip boundary condition on the heliopause at large radii. The rotational velocity about 
the heliotail in the near-field solution decays as the fourth power of distance from the shock. The polar 
divergence of the streamlines will have consequences for the previously discussed magnetic pressure ridge 
that may build-up just inside the heliopause. 

I. INTRODUCTION 

We have previously shown the conditions under which potential flow occurs beyond the terminal shock of the 
solar wind. Topologically, potential flow represents well many of the global properties which are universal to 
all heliospheric stagnation point flow topologies [5 m ess, 1990]. The flow is self-consistently incompressible to 
about 15% and irrotational to 10% [Suess and Nerney, 1990, 1991 (SN1, SN2)]. The latter is true when the 
non-slip condition is used on the heliopause at large radii. Because it is likely that reconnection occurs near 
the stagnation point between the heliosheath and interstellar magnetic fields, we believe that the heliopause 
is best represented as a non-slip interface [Suess and Nerney , 1993], This reduces the solar wind Mach 
number down the heliotail compared with free-slip numerical models, reducing the effect of compressibility 
in this region. 

Recently, a numerical simulation for an external Mach number of M=0.8 has shown that the shock is nearly 
spherical and Sun-centered even in that extreme case [ Steinolfson , 1994]. That calculation further showed 
that the density changes by a total of only 28% in crossing the heliosheath and heliopause and that, in 
general, compressibility is not globally important. The only place where it is not completely negligible is in 
the boundary layer at the heliopause. Vorticity is introduced only in a thin layer just beyond the terminal 
shock. This lends credibility to potential flow modeling in the case of a subsonic external flow. We do not 
consider the case of external supersonic flow because we believe that the interstellar wind is likely to be 
sub-Alfvenic and information will be able to propagate upstream, removing the necessity for a bow shock 
(also see the discussion in Nerney et a/., 1994). 

The merging of the shocked, subsonic solar wind with the 25 km/s interstellar wind was calculated in our 
potential models by including the effects of a finite radius ( R s ) terminal shock. This effect entered the 
solution for the streamlines through a parameter e which is a measure of the ratio of the interstellar ram 
pressure to that of the pressure in the shocked solar wind just beyond the terminal shock. Vanishingly 
small values of e reproduce Parker's [1963] model for a weak interstellar medium with a large turning radius 



for the solar wind. Because of this, the shock radius is small compared to the transverse dimension of 
the heliosphere, and the shock appears to have a negligible radius. The numerical solutions were plotted 
for variable e, and the distance from the stagnation point to the terminal shock was no larger than R s for 
probable heliospheric parameters. A plot of the streamlines illustrating this result for e=.125, an appropriate 
value for the heliosphere, is shown in Figure 1, which also serves to define the related terminology. This 
problem is of general astrophysical interest for the merging of shocked stellar winds with the interstellar 
medium. 

It is our purpose now to show the nature of the change in these results when a heliocentric latitudinal ram 
pressure variation is imposed on the supersonic wind upstream of the shock. This requires a non-spherical 
shape for the shock due to, for instance, high-speed streams emanating from the poles of the sun, and changes 
the character of the solutions near the shock. Observations of a latitudinal variation in the interplanetary 
Lyman a distribution [Kumar and Broadfoot, 1979; Lallement et a/., 1985] have been interpreted as indicating 
a 30-50% decrease in the solar wind mass flux from solar equator to pole during 1973-1977. This would lead 
to a increase in solar wind ram pressure from equator to pole in the outer solar system for a 30% decrease in 
mass flux, but not for a 50% decrease. Our model is more appropriate for periods with mass flux variations 
in the lower range of values. 

We have linearized the solutions in the departure from sphericity. The mapping from heliocentric spherical 
polar coordinates to heliocentric heliotail coordinates (cylindrical) introduces an azimuthal dependence of 
the streamlines about the heliotail. There is also rotational flow about the heliotail which decays as the 
fourth power of distance from the shock. A further asymmetry is introduced because the footpoints for each 
streamline begin at different distances from the sun. 

We present the mathematical formalism in II, the numerical analysis in III, and the conclusions in IV. 

II. MATHEMATICAL FORMALISM 
The potential flow equations are the following: 


V * v = 0 

(la) 

v — — V<F 

(16) 

V 2< J> = 0. 

(lc) 


The boundary conditions are: 



e r • v{9,<j>) \ R , = v 8 (0,<j)) (2a) 
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where R s is the distance to the terminal shock, v 3 is the solar wind flow velocity downstream of the shock, 
<j> is the azimuthal angle about the heliotail, and 0 is the polar angle with respect to the z-axis, the direction 
of the interstellar wind. We impose the following asymptotic flow pattern on the supersonic solar wind: 


Vs = v 8 o(l - 6 sin 2 0 5 ) (3a) 

where 6s is the solar polar angle and 6 < .2, and t; 5 o is the solar wind velocity downstream of the shock 
along the solar rotation axis. Rotating into heliotail, heliocentric spherical polar coordinates: 


Vs = v s0 [l — 6(cos 2 6 -f sin 2 0cos 2 <£)], (36) 

The rotation of coordinate systems introduces a cos 2 <j) dependence on azimuth about the heliotail. Equiva- 
lently, 


v $ = v,o[l - ~(4P 0 + 2P 2 + cos 2<f>)] (3c) 

where Pi and P™ are the Legendre polynomials and the associated Legendre polynomials, respectively. 

Now, <& may be written as: 


*=£ 


/ r \ n / T> \ n + l 

Am,n \Ro J cos m<j> P™ (cos 0) 


(4) 


and R s o is the value of R s along the solar rotation axis. The A m>n are determined by requiring that the 
velocity of the shocked wind approach that of the interstellar medium at large r, eq.(2b). As in SN1 


A m>n — 0 V (m, n) 


(5a) 


except that 


A),l — ^ioo^iO 


(56) 


where 


Rs 


R»ov,{6,<l>) _ R ^ v 


( 6 ) 


which is derived by setting the ram pressure of the solar wind equal to the net pressure in the local interstellar 
medium (LISM). Eq. (6) serves to define V and gives the angular variation for the distance to the termination 
shock. 



Taking the radial derivative of equation (4), including only terms present in v ,, and evaluating at R a , 
find: 

^ C I 3 , r»0 , p r>2 
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’ cos 2<f>) 


we 


( 7 ) 
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so that the angular variables occur in both the numerators and denominators of equation (7), unlike the 
spherical shock solutions. We have defined e as in SN2: 


e — 


V*0 


3/2 


(8) 


Now equating equations (3c) and (7) and linearizing the various powers of l/V using binomial expansions, 
allows us to solve for the coefficients of the spherical harmonics so that $ may be written as: 


$ 
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where 


and 


r 2 = w 2 + z 2 


P 2 = o (3 COS 2 0—1) 
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Pj = 3 sin 2 I 


The velocities are derived from the derivatives of 4>. 
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In passing, we note, as expected, that there are no order e effects in 


We will now switch to cylindrical coordinates aligned with the heliotail and convert the derivatives in (11) 
to: 
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For economy of notation we have written r 2 for w 2 + z 2 . The streamlines are defined by 

dz dzu zod<j> 


v z v. 
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Using the nondimensional variables from SN1 
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(note that € X ^ 3 /R s q remains finite as e — ► 0) the streamline equations can be written as: 
dC drj 


d<j> 
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where 
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and Fi , F 2 are rewritten as 
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(14a) 
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When rj = 0 and drj/d( is infinite, v z must be zero. Setting the left denominator of (15) to zero then 
determines the non-dimensional distance to the stagnation point: 


<c(l - 6F 2 ) + (C 2 ) 3/2 - c( 1 - 2.46) = 0 


(19) 


III. NUMERICAL ANALYSIS 


Figure 2 shows the streamlines for a 20% deviation from a spherical shock both in the meridional plane (solid 
lines; 0 — 90°, which includes the solar rotation axis) and in the equatorial plane (dashed lines; 0 = 0°). 



The equatorial streamlines map back to the inner circle, a cross-section of the terminal shock, while the 
meridional streamlines map back to a cross-section of a prolate spheroid and show the 20% bulge over the 
poles. The positions of the heliopause for the polar and equatorial streamlines are indicated on the right 
edge of the graph, showing that the bulge in the shock over the poles pushes the heliopause out in the 
meridional plane compared to the equatorial plane. The streamlines are nearly radial on the outside of the 
shock for <t> = 90°, but there is elevated non-radial flow on the shock for <j> = 0° and e = .125. We choose 
to superimpose the graphs of the streamlines for <j> = 0°,90° because these are the two angles about the 
heliotail for which v<f> is zero (see eq. (lid)). The following plots indicate that the problem is inherently 
three-dimensional for any other value of <j>. 

Figures 3a, b, and c show three views of streamlines that begin on the terminal shock at <j> = 60°, followed 
by fig 3d which superimposes both <j> = 30° (grey curves) and <j> = 60° (black curves). The effects of v<f> 
are seen in the near-field solutions for the potential flow but decay asymptotically as r -4 (see eq. (lid)). 
The potential flow solutions connect the axisymmetric boundary condition at infinity (all streamlines blend 
into the assumed uniform interstellar wind) with the anisotropic boundary condition on the terminal shock 
(high-speed streams at the poles) so that the streamlines must diverge away from the poles. Another way to 
describe this effect is to remember that the solution to Laplace’s equation minimizes the total kinetic energy 
per gram integrated over any volume bounded by streamlines. This requires a divergence of the flow away 
from the poles where the kinetic energy per gram on the shock is high. 

IV. SUMMARY AND CONCLUSIONS 

We have solved a potential flow model for the flow down stream from a non-spherical terminal shock. The 
equations were linearized for small departures from sphericity and the streamline equations were numerically 
integrated to show the three-dimensional flow effects. The anisotropic streamlines on the terminal shock 
diverge away from the polar regions, while in the distant heliotail they smoothly take on the characteristics 
of the assumed uniform streamlines in the interstellar medium. This would have consequences for the 
previously reported ridge of magnetic flux that may build-up just inside of the heliopause [Nerney, Suess , 
and Schmahl , 1993]. The divergence of the flow away from the poles will initially open magnetic flux tubes, 
reducing the magnetic flux in the ridge. We have also reported the magnetic field line topology and solar 
cycle imprint for a kinematic magnetic field in the heliosphere [Nerney, Suess , and Schmahl , 1994]. A non- 
spherical shock would change the detailed shapes of both the magnetic field lines and the solar cycle imprint 
shown in that calculation. However, the largest difference shows up in using a free-slip versus a non-slip 
boundary condition on the heliopause. In the present calculation we have applied a non-slip boundary 
condition. The consequence of a free-slip boundary condition is to increase the flow speed down the heliotail, 
with a corresponding decrease in the cross-sectional area of the tail and in the thickness of the heliosheath 



on the upstream side of the heliosphere. 
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FIGURE CAPTIONS 


Figure 1. Streamline plot for the potential flow solution [Suess and Nerney, 1990], for e = .125 - an 
appropriate choice for the heliosphere. The termination shock, stagnation point, heliopause, heliosheath, and 
heliotail are indicated. The stagnation point of the flow is at R c = 2 R s for this value of e. The cylindrical 
coordinates are labeled for distance down the heliotail, z, and cylindrical radius about the heliotail, w. The 
flow is symmetric about the heliotail (independent of <j>) for a spherically symmetric shock, so the streamlines 
are the same in any plane containing the axis of the heliotail. 

Figure 2. A plot of streamlines for the equatorial plane (<f> = 0°; dashed lines and circular shock) and the 
meridional plane ( <j> = 90°; solid lines and oval shock for a 20% elongation over the poles). Arrows on the 
right-side of the figure (A for <j> — 0°, Bfor90°) point to the respective heliopauses. Streamlines are different 
in other planes of constant <j>. 

Figure 3. Panels (a), (b), and (c) show different views of the streamlines that begin on the terminal shock at 
<f> = ® 0°, while (d) shows a superposition of <j> = 30° streamlines (grey curves) and <f> = 60° streamlines (black 
curves). Only panels (b) and (c) include streamlines in the interstellar medium. In panel (a) the view is 
looking directly at the sub-stagnation point on the terminal shock, while in (b) and (c) the stagnation point 
is to the right. The streamlines diverge from the poles as the anisotropic boundary condition on the shock 
(high-speed streams at the poles) evolve into the uniform flow of the interstellar wind in the downstream 
heliotail. 
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